Contact forces in regular 3D granular pile 
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Abstract 

We present exact results for the contact forces in a three dimensional static piling 
of identical, stiff and frictionless spheres. The pile studied is a pyramid of equilateral 
triangular base ("stack of cannonballs" ) with a FCC (face centered cubic) structure. 
We show in particular that, as for the two dimensional case, the pressure on the base 
of such a pile is uniform. 

PAC numbers: 

83.70.Fn: Granular solids 

46.10.+z: Mechanics of discrete systems 

The force propagation through a heap of granular materials is a controversial issue (see 
refs. [|], and references herein) and still attracts much attention. Experiments in conical 
sand piles show the existence of a depression in the normal force on the base (a dip) under 
the pile's apex rather than the intuitively expected maximum |3], |J. Much work has been 
done in this subject in order to explain this effect numerically and theoretically . Most of 
the results obtained in this domain are for models considering two dimensional (2D) piles 
of discs or rods || ||, 0, ||, ||, [l(| [H], |12| . It is unclear however, whether in general, these 



results can be directly generalized to three dimensions (3D). In the case of the continuum 
approach proposed in ref. [|F| it was noticed that the choice of the second closure equation 
needed for the 3D case did not have an important impact on the results obtained, so that 
the 2D results survive the passage from 2 to 3 dimensions. This case is a priori unlikely 
when discrete models are considered. One difficulty might be that a simple 3D extension 
of a 2D pile, which is obtained by a revolution around the symmetry axis of a 2D pile 
is not any regular 3D piling of spheres but rather a 3D piling of tori. Another possible 
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problem is that in 3D there are more packing possibilities than in 2D. Even if one accepts 
that each layer of spheres is arranged on a triangular lattice, there are different ways to 
pile layers one on top of the other, giving different arrangements; FCC, HCP, DHCP etc. 
Another issue is the shape of the base. A circular base, as for the real conical sandpile, is 
"unnatural" for spheres in the 3D context, while in 2D it is the simplest one since for a 
2D pile, a circular base is simply a segment. 

Here we study a 3D case that might be considered as an extension of the one proposed 
in refs. \TA, 12] for 2D. The configuration we study is presented in fig. |l|. It is a regular 3D 



pyramid with triangular base and an FCC-like structure^, composed of identical, spherical, 
stiff and frictionless particles of mass m and radii R = 1 (this piling is also known as 
the "spanish stack of cannonballs" ) ||15|| . In each horizontal layer the distance between 
neighboring spheres is 2a, where l<a<3/2so that there is no overlap between spheres, 
horizontally or vertically. Our aim is to calculate the contact forces everywhere in the pile, 
especially the normal and shear forces applied by the pile on the supporting surface. For 
the reduced 2D case it was shown |L4| that the pressure profile on the supporting surface is 



uniform. This feature was also proved in ref. []TB| for a similar pile with periodic vacancies. 

In the following, we will index the spheres using a triplet of integers k) where k is 
the horizontal layer index counted from the top of the pile down (the topmost layer has 
k = 0) and the couple indices the spheres in the k-th layer as shown in fig. |2|. If the 
pile has a total of L layers then, < k < L— 1, < j < k and < i < k — j. Components 
of the force vectors will be given in the vectorial base shown in fig. [I]. 

Let us take a closer look at one of the spheres, indexed: (i, j, k). We denote, as shown 
in fig. [5], by a, b and c the forces applied by the downwards neighboring spheres, and by 
a', b' and c' the forces applied by the upward spheres. The weight of the particle is w. 
One can easily conclude from fig. ^| that 

a = ae a , a' = — a'e a with e a = (sin#sin0, sin#cos0, — cos#), (la) 
b = /3e&, h' = —(3'eb with e& = (sin # sin 0, — sin # cos 0, — cos6>), (lb) 
c = 7e c , c' = — 7'e c with e c = (— sm.9, 0, — cos#), (lc) 



where a, ft, 7, a', ft', 7' are the norms of the vectors a, b, c, a', b', c' respectively, and e a , eb, e c 
are the unit vectors along the the vectors a, b, c respectively. 

Since the pile is on top of an horizontal surface the weight of each sphere is give by 

W = 77^(0, 0,1). (2) 

The equilibrium of each sphere is given by 

^F = a + b + c + w = 0. (3) 

1 the structure is exactly FCC if a = R, in which case the diagonal of the elementary cube is oriented 
vertically 
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Combining eqs. and [3] and simplifying we get the following linear system: 



a sin — a sin + (3 sin — (3' sin — 7 + 7' = (4a) 

a-a'-p + fi = (4b) 
a cos 9 — a' cos + /3 cos 6 — (3' cos 6 1 + 7 cos # — 7' cos = mg (4c) 

Giving finally 

Aa + A/3 - A7/sin0 = (5a) 

Aa-A/3 = (5b) 

Aa + A(3 + A>y = mg/ cos 0, (5c) 

where Aa = a — a', A/3 = (3 — [3' and A7 = 7 — 7'. This system is easily solved yielding 

a a n my 1 a sin0 mo . . 

Aa = Ap= — and A7 = — - — (6) 

2(1 + sin 0) cos 1 + sin cosy 

In the case of an equilateral triangular base we have = 7r/6, so that 

Aa = A/3 = A 7 = (7) 
3 cost' 

In order to get the force acting on a specific contact one should count the number of 
spheres up to the free surface in the corresponding direction, each one of these spheres 
contributing mg/(3cos9). We denote N a (i,j,k), Np(i,j,k) and N y (i,j,k) the number of 
particles on top of the k) sphere in the corresponding direction, so that 



a'(i,j,k) = mg/ (3 cos 9) N a (i, j,k), (8a) 

P\i,j,k) = mg/(3cos9) N p {i,j,k), (8b) 

iXhhk) = mg/ (3 cos 9) N 7 (i, j, k) . (8c) 

One can easily show that (see appendix for details) 

N a (i,j, k) = k-i-j, Np(i,j, k) = i and N y (i,j, k) = j, (9) 

yielding, when introduced into eq. |c| 

a'(i,j, k) = (k — i — j)mg/{3 cos 9) (10a) 

(3\i,j,k) = img/ (3 cos 9) (10b) 

l'(i,3,k) = jmg/ (3 cos 9), (10c) 
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which gives us the amplitude of the contact forces anywhere in the pile. If we suppose that 
the pile is L layers high and that the bottom layer spheres are "glued" to a flat surface, the 
force applied by the bottom layer spheres on the surface; F s (i,j) can easily be calculated, 
since from the equilibrium of a bottom layer particle we get 

F a (i,j) = a!(i,j, L - 1) + b'(i,j, L - 1) + c'(i,j, L - 1) + w, (11) 

and using eqs. [L],|2| and pi] in this last equation gives 

F s {i,j) = f -1/2 (£-1-3.7) sin^,- V3/2 (L - 1 - 2t - j) sin 9, 

3 COS 9 \ 

(L + 2) cose). (12) 



Thus, the normal force on the base of the pile applied by each one of the spheres is the 
same and has an amplitude of 

mg(L + 2)/3. (13) 

This force is independent of 9 and hence, of a. One can easily verify that the sum of the 
normal forces on the surface is equal to the total weight of the pile mg L(L + 1)(L + 2)/6. 
The shear force applied on the surface is given by 



— mg tan 9 



((L - 1 - 3j)/2, V3/2 (L — 1 — 2i — j)) . (14) 



An example of this vector field is given in fig. (|. Notice that the shear force vanishes under 
the apex of the pile as required from symmetry. 

In conclusion, in the case of a FCC-like piling with equilateral, triangular base of stiff, 
identical, spherical particles we were able to calculate the contact forces between any couple 
of neighboring particles (eq. [K|). In particular we have shown that normal forces on the 



base are constant and are given by eq. [13] while the shear force is given by eq. [14. These 



results are similar to those obtained in ref. |14| for the 2D case and might suggest some 



applicability of 2D model results to 3D or, at least, that the starting point in 3D is the 
same as for 2D. It is clear, however, that horizontal contacts should be considered in order 
to get more physical behavior like arching. But in this case, the system is hyperstatic and 
more complicated techniques are to be used, exactly |], Q or numerically J7|, in order to 
get the correct contact network. Further investigations should be carried out in order to 
obtain similar results for other base shapes like hexagons and other structures, especially 
in the RHCP (Random Hexagonal Compact Packing) case that can be a convenient way 
to introduce disorder in the system so to approach a modelization of real granular piles. 
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Figure 1: A pile of four layers with a = 1.2 R. 
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Figure 2: Indexing convention of the spheres for a layer in the pile. The layer shown is 
the one for k = 3 (forth from the top). Remark that force components are given in the 
vectorial base shown in fig. |l| and not in the one used for indexing. 
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Figure 3: A schematic drawing of the forces acting on a given sphere centered at O (not 
shown). The neighboring spheres are the ones centered at A, B and C (layer beneath) 
which apply the contact forces a, b and c respectively, and A',B' and C (layer above), 
which apply the contact forces a', b' and c' respectively, w is the weight of the particle. 
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Figure 4: A plot of the shear force applied by the pile on the supporting surface for a 10 
layer pile. Notice that the shear force vanishes under the apex of the pile as required by 
the symmetry of the pile. 
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1 Appendix 



Let us start by showing by induction over k that : N 7 (i,j, k) = j V{0 <j<A:|0<2< 
k — j}. When k = it is clear that iV 7 (0, 0, 0) = since the top particle does not have 
any neighbor on it. Now let us suppose that N 7 (i,j, k) = j V{0 < j < k\0 < i < k — j} 
for a certain k and prove it for k + 1. If j = the particle is on the surface of the pile and 
iV 7 (i, 0, k + 1) = = j. When j > the (i, j, k + 1) sphere is in contact with the (i, j — 1, k) 
particle, so if k + 1) particle is centered at point O in fig. |], the {i, j — 1, k) particle 
is centered at C and the contact point is at the middle of the segment OC. Thus, we 
conclude, using the hypothesis, that N^(i,j, k + 1) = N 7 (i,j — 1, k) + 1 = j, Q.E.D. 

N a (i, j, k) and Np(i, j, k) are obtained by a rotation of the pile around the z axis of 2ti/3 
and 47r/3 respectively, which correspond to applying the transformation {j — > k — i—j, i — > 
j} once and twice respectively. 
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